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Abstract 

In the setting of the principle of local equilibrium which asserts that 
the temperature is a function of the energy levels of the system, we ex¬ 
hibit plenty of steady states describing the condensation of free Bosons 
which are not in thermal equilibrium. The surprising facts are that the 
condensation can occur both in dimension less than 3 in configuration 
space, and even in excited energy levels. The investigation relative to 
non equilibrium suggests a new approach to the condensation, which al¬ 
lows an unified analysis involving also the condensation of g-particles, 

— 1 < q < 1, where q = ±1 corresponds to the Bose/Fermi alternative. 
For such q-particles, the condensation can occur only if 0 < q < 1, the 
case 1 corresponding to the standard Bose-Einstein condensation. In this 
more general approach, completely new and unexpected states exhibiting 
condensation phenomena naturally occur also in the usual situation of 
equilibrium thermodynamics. The new approach proposed in the present 
paper for the situation of 2 nd quantisation of free particles, is naturally 
based on the theory of the Distributions, which might hopefully be ex¬ 
tended to more general cases. 
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1 introduction 


The investigation of the Bose-Einstein Condensation (BEC for short) has a very 
long history after discovering a new particle statistics by Satyendra Nath Bose 
(cf. [10]) and Albert Einstein (see e.g. [14]) at the beginning of 20th century. 
It concerns the fact that a macroscopic amount of elementary particles having 
integer spin (called Bosons in honour of N. Bose) can occupy the ground state 
after the thermodynamic limit. The simplest and most important ideal model 
is that consisting of free massive Bosons. We also mention the Gross-Pitaevskii 
equation (cf. [21, 35]) describing a weakly interacting Bose-Einstein condensate 
at low temperature. Many phenomena involving quasi-particles like phonons 
or magnons can be described by the BEC. Recently, in [25] the condensation 
of massless particles like photons has been pointed out as well. Also the phe¬ 
nomenon of the superfluidity of the helium isotope He 4 seems to be tightly but 
not directly connected with the BEC. In fact, on one hand the superfluidity is 
mentioned as one of the most important examples confirming the evidence of 
the condensation effects. On the other hand, it is unclear if this phenomenon is 
directly connected with the condensation of (quasi) particles, and was explained 
by Lev Landau (cf. [27]) with a fluid composed of two indivisible components, 
one in the superfluid phase and the other one in the normal phase. At zero 
temperature, the He 4 is entirely composed by the superfluid component. The 
thermal agitation causes the excitation of the part of the dispersion spectrum 
corresponding to phonons, and mainly to rotons. In such a way, as the tem¬ 
perature increases, the portion of the superfluid component becomes lover and 
lover until to disappear at the critical temperature, known as the A-point. 

The behaviour at low temperatures of particles obeying to Fermi-Dirac 
statistics is completely different. Due to Pauli exclusion principle, it is well 
known that Fermions (i.e. quantum particles of half-integer spin) does not lead 
to any condensation by the Pauli exclusion principle. Nevertheless, the other 
isotope helium He 3 still exhibits superfluidity at a temperature very close to 
0° Kelvin, even if these are Fermi particles. The superfluidity of He 3 can be 
justified with the fact that at very low temperature, H e3-particles form the so 
called sea of pair-particles which can be considered as Bosons. Namely, also 
He 3 can exhibit superfluidity. According to the BCS theory (cf. [7]), it is 
precisely the same phenomenon occurring in superconductors where pairs of 
electrons forms the so called Bardeen-Cooper pairs. The reader is also referred 
to [18, 19, 20], where it is shown the condensation of Bardeen-Cooper pairs in 
arrays of Josephson junctions. In this sense, the phenomenon of superconduc¬ 
tivity is also connected with the BEC. Finally, in Sections 11.7 and 11.10 of the 
monograph [6], the BCS pairing appears as a general phenomenon of Bosonisa- 
tion in the setting of stochastic limit. For a exhaustive treatment of the various 
condensation phenomena and perspectives, we mention the sample very far to 
be complete, of books [29, 32, 37], and [13] with the literature cited therein, for 
a rigorous mathematical approach to BEC. 

Typically BEC is considered in states which are in equilibrium at a certain 
temperature and with respect to the natural dynamics generated by a given 
Hamiltonian. The definition of equilibrium states in terms of Kubo-Martin- 
Schwinger (KMS for short) condition (cf. [26, 33]) w.r.t. to a given dynamics 
has the advantage compared with the usual Gibbs prescription given in terms 
of the Hamiltonian implementing the dynamics, of being valid also after the 
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infinite volume limit, see e.g. [22]. It is known that an equilibrium state for a 
given dynamics is automatically stationary for this dynamics. 

Non Equilibrium Steady States (NESS for short) play also an important role 
in non equilibrium thermodynamics. A natural class of NESS naturally emerges 
from the stochastic limit of quantum theory, see e.g. [2, 3, 5, 6] and references 
therein. States in this class are called local equilibrium, (or local KMS states) 
where the term local refers to the fact that the temperature is a function of 
the Hamiltonian of the system so that any energy level can be thought to be 
in equilibrium at a temperature depending on the level itself. Moreover, these 
NESS are characterised in terms of a Local Equilibrium Principle, that can be 
considered as a generalisation of the usual KMS equilibrium one. The Local 
Equilibrium Principle was formulated in [2] for systems with pure-point spec¬ 
trum dynamics. In [1] this principle was formulated for systems whose dynamics 
is implemented by a strongly continuous 1-paranreter unitary group acting on 
a given Hilbert space. Several equivalent (or almost equivalent) formulations 
of it were also discussed for systems with pure-point spectrum dynamics. The 
present paper extends the above results in several directions as follows. 

(i) Exploiting the fact that the quasi-free dynamics (i.e. one-parameter group 
of Bogoliubov automorphisms) leave the Weyl C*-algebra invariant, we give a 
purely LA-algebraic formulation of the Local Equilibrium Principle for quasi- 
free (gaussian) NESS with respect to a given quasi-free dynamics, and a given 
inverse temperature function 0. 

(ii) We characterise the structure of these states and we find that, for the affine 
function 0(h) = 0(h)h = '-fjf- (with the Boltzmann constant k = 1), hence 
including the chemical potential p, they are reduced to the usual quasi-free 
equilibrium states. 

(iii) We extend the study the BEC phenomenon to the setting of local equi¬ 
librium quasi-free states that, because of (ii), includes the equilibrium ones as 
particular cases. 

(iv) In the setting of local equilibrium, we exhibit states exhibiting BEC, not 
only on the ground state, but also (or only) on some excited levels, depending 
on the function 0, entering in the definition of the Local Equilibrium Principle. 

(v) We construct examples of local NESS for which BEC occurs also for spatial 
dimensions d different from the usual ones (d > 3 for free massive Bosons, and 
d > 2 for massless particles). 

(vi) We construct local quasi-free NESS, exhibiting BEC in excited levels, for 
which the rotation symmetry can be spontaneously broken. The rotationally 
invariant ones are obtained by averaging on the spheres in momentum space 
[p € R d 0(h( p)) = 0}, whenever h( p) is the one-particle Hamiltonian of a 
free massive Boson. 

(vii) We extend our analysis to include the case of q-Defornred Commutation 
Relations with q £ [—1,1], and we prove that BEC can occur also for the q- 
relations provided that q € (0,1]. The case q £ (0,1) opens some mathematically 
(and maybe also physically) interesting possibilities which at the moment seem 
to be unexplored. 

(viii) Even in the usual case of equilibrium thermodynamics, we exhibit com¬ 
pletely new states describing BEC, which are mathematically meaningful, and 
could have promising physical applications. 
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Many of the results of the present paper concerning the 2 nd quantisation of 
free particles, are obtained by using a new approach based on the theory of the 
Distributions, which hopefully might be extended to more general situations. 


2 the local equilibrium condition 

As a first step in posing the problem, we consider a physical system with a finite 
degrees of freedom, whose observables are typically described by the C *-algebra 
of all the bounded operators 23 (2f) acting on a separable Hilbert space 2C. The 
time-evolution is given in Heisenberg picture by 

a G ®(5f) ^ at. (a) = e lHt ae~ lHt G ®(2f), <Gt, 

where the Hamiltonian H > 0 of the system is a self-adjoint operator acting on 
2C, which is generally unbounded. To simplify the matter, we suppose that 

OO 

H = '^£ n \i’n){i’n\ ( 2 . 1 ) 

71=0 

is positive with compact resolvent, with eigenvalues e n G a{H) and the corre¬ 
sponding eigenvectors tf> n G Jf, repeated according the multiplicity and rear¬ 
ranged in increasing order if needed. Fix any positive function (i.e. the local 
inverse temperature ) (3 : R+ —» R+ such that is a trace-class operator. 

According to the Local Equilibrium Principle (see e.g. [1] for details and the 
literature on the topic), define the state 

ujp{a) := Tr (e~^ H)H a) , a G ®(5C). (2.2) 

Denote 3 r (51) the sub algebra consisting of all the finite-rank operators acting 
on 2f, and Vij := \ipi)(ipj\, the i/’j being the eigenfunctions of H according to 
(2.1). It is immediate to show that 

zGCg uip(aa z (b )), a G 23 (J{), b G T(2t), 

is well defined and entire. The Local Equilibrium Principle (LEQ for short) 
simply means 

up(aa t +ip( Ei ){vij)) = e [/3(ei)_/J(ej ' )l£j w / 3 (at(u ii )a), a G $(21). 

The LEP, which reduces to the usual KMS boundary condition (cf. [13, 22]) 
when [3 is the constant function, is not immediately generalizable to arbitrary 
dynamical systems (21, at) consisting of a C*-algebra, and an action a t of one- 
parameter group of possibly outer ^-automorphisms of 21. A natural way to get 
such a possible generalisation at least for inner time-evolution, is to look at the 
modified evolution 

a G 'B(Jf) ^ af\a) := e *P(H)Ht ae -zP(H)Ht e (2.3) 

In the discrete spectrum Hamiltonian case discussed below, 

z€C^Ufi(aa^(b)) f a G 23(J{), b G S^Jf), 
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is again well defined and entire, and the state top satisfies the KMS boundary 
condition w.r.t. this modified evolution at inverse temperature 1/T = 1. 

We end by reporting the property analogous to the well-known one concern¬ 
ing the KMS boundary condition in momentum space. Consider any continuous 
compactly supported function / on K, together with its Fourier (inverse) trans¬ 
form /. Then for each a £ 'B(J-C) 

/ +oo c+oo 

f(t)up(a*a ( /\a)), i/ a {f):= f(t)up(a[ 0) (a)a*) 

-oo J — oo 

define Radon measures on R. In this simple situation, /i and v are linear com¬ 
binations of Dirac measures supported in a subset 


depending of a. In addition, such measures are equivalent with Radon-Nikodym 
derivative given by 


d^q 

d^q 


(k) = e~ k 


where the k £ K a . The last simple calculation has the following meaning. 
The set the Bohr Frequencies of a physical system as before consists of all the 
possibly energies of the quanta that it can emit or absorb while interacts with the 
environment, see [4] for a simple explanation. It is given by {s n — e m \ e n , £ m £ 
a(H)}. In the more general case of LEQ, the natural Bohr Frequencies (or 
the so-called Arveson Spectrum in the mathematical language) for the possible 
transition of the system are precisely those {/3(e„)e n —/3(e m )£ m | e ra ,£ m €E <r (H)} 
of the modified Hamiltonian Hp := (3(H)H. We end the present section with 
the following facts. First, we observe the similar dependence of the temperature 
directly on the Arveson Spectrum of the involved automorphism arising from 
the spectrally passive states in [15]. Second, we point out again quasi-free states 
with non constant temperature described in [12] for non equilibrium relativistic 
thermodynamics. 


3 local equilibrium principle for free Bosons 

Let f) the one-particle Hilbert space. The Canonical Commutation Relations 
(CCR for short) algebra describes Bose particles and is generated by 

a(f)aHd) ~ a\g)a(f) = (f\g)l f,g£ f) 0 . (3.1) 

for a dense subspace f)o C f). It is well-known that the CCR cannot be repre¬ 
sented as bounded operators acting on some Hilbert space even if f) = C. The 
standard procedure is to pass to the CCR represented in the Weyl form. It is 
the universal C*-algebra CCR(bo) generated by unitary operators satisfying the 
relations 


W(h)* = W(—h) , W(0) = 1 , (3.2) 

W(f)W(g)=e ilJ ^W(f + g), f,g£bo. 

Thus, with CCR(f)o) we denote the universal C*-algebra generated by the Weyl 
operators {W(f) \ f £ f)o} satisfying the relations (3.2). 
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Let h > 0 be the one-particle Hamiltonian. Suppose that e lht bo = f)o- 
Then the time evolution is generated by a one-parameter group of Bogoliubov 
automorphisms which acts on the Weyl generators as 

a t (W(f)) :=e ztdrw W{u)e- itdrw = W(e lth f), /Gf) 0 ,tGR. (3.3) 

Here, H = dT(h) is the second quantised Hamiltonian, and we have identified 
W(f) with its image 7rp(PT(/)) in the Fock representation because CCR(t)o) is 
a simple C*-algebra, see e.g. [13, 30] and the references cited therein. 

For systems for which 1) is finite-dimensional, the second quantised Hamil¬ 
tonian H := dL(ft.) has compact resolvent, which typically occurs when one 
considers a free gas confined in a box of finite volume. In this framework, 
e -p(H)H j g automatically trace-class, see e.g. Proposition 5.2.27 of [13]. Thus, 
the unique state uip satisfying the LEP w.r.t. the time evolution (3.3) is the 
quasi-free state uniquely determined by the two-point function 

u p (a\g)a(f)) = (f\(e^ h -l)- 1 \g) i . f,g€ ft. 

By using the CCR (3.1), we have 

W /3 {a(g)a i (e-^ h)h f)) = w p (a^{f)a{g)) , (3.4) 

{g)a(e 0{h)h f)) = wp(a(f)a f (g)). 

The boundary conditions (3.4) gives the natural generalisation for LEP in the 
cases of quasi-free states of CCR algebras, providing several nontrivial examples 
for which the dynamics is neither inner, nor has pure-point spectrum in general. 

We start with a separable Hilbert space f), together with a self-adjoint 
positive Hamiltonian (in general unbounded) h > 0, which are nothing but 
the one-particle Hilbert space with the relative one-particle Hamiltonian. Let 
EGlq e(e) £ 23 (J~C) be the right-continuous resolution of the identity associ¬ 
ated to h (cf. Section VIII.3 of [36]) according to 

h = Jedehie)- (3.5) 

Fix a Borel function /3 : M + R + which is positive almost everywhere w.r.t. 
the measure class determined by the projection-valued measure de(e). The key- 
point to manage open thermodynamical systems where exchange of matter is 
also allowed, is to introduce the chemical potential, see e.g. [28]. In order to 
achieve the chemical potential for local equilibrium, we follow the fact pointed 
out in Section 2 which asserts that a state satisfying the LEP is indeed an 
equilibrium one for a modified Hamiltonian at inverse temperature 1/T = /3 = 1. 
For the activity z (cf. pag. 47 of [13]), it leads to z(l/T, p) = z( 1, p) = e M , where 
p is the chemical potential. This justifies the definition for p £ K. in 

7/i(*) “ e /3(x)x -' 1 • (3.6) 

For the g-particles with q £ (0,1] we have a natural restriction to the value of 
the allowed chemical potentials p coming from the positivity condition of the 
occupation number, see below. For the Bose case q = 1, such a condition simply 
leads to p < 0. 
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A quasi-free state is a mean-zero gaussian state to, a-priori defined on the 
whole CCR(f)). Is uniquely determined by its two-point function as it is ex¬ 
plained in Section 5.2.3 of [13]. In fact, we can start with a quadratic form Q 
on f) such that 

WQ(a f (/)a(/)) =Q{f), /el). 

This means nothing but that in a more rigorous form, 

u Q (W(f)) = exp{—(||/f/4 + Q(f)/ 2)} , / 6 h , 

which is meaningful also when Q(f) = + 00 . A standard procedure to manage 
quasi-free states (see Theorem 5.2.31 of [13], for the situation describing BEC) 
is to consider the domain Dq C 1) made of all the elements on which Q is finite, 
and consider the sesquilinear form Fq defined by polarisation on all of Dq such 
that 

UQ(a\f)a{g)) = F(gJ) , /, <? € Dq . 

Thus, we fix a value of g < 0, and a dense subspace l)o C I) such that 

(i) e lh % = l)o, t € R. 

Among the quasi-free states ujq on CCR(l)o), we consider those for which 

(ii) f)o C Dq, and 7 /i (/i)f)o C Dq, with 7 ^ in (3.6). 

The natural generalisation of the LEP to quasi-free states of Bosonic systems is 
one of the (equivalent) conditions given in (3.4). 

Definition 3.1. Suppose that the dense subspace f)o C I) and the quasi-free 
state u> € S(CCR(f)o)) fulfill (i)-(ii) above. Then w satisfies the Local Equi¬ 
librium Principle at local inverse temperature f3 and chemical potential g < 0 
if 

u{a\f)a{^(h)g)) = u>(a(g)a\f )), f,g € f ) 0 . (3.7) 

It is immediate to prove (cf. Section 5.3.1 of [13]) that a quasi-free state 
satisfying the LEP is invariant w.r.t. the modified evolution (2.3), provided 
etPWhtfo = fj 0) t e R. In general, it is not clear if it is invariant also w.r.t. 
the natural time-evolution of the system generated by the one-parameter Bogoli- 
ubov automorphisms 1 1 —► e lht . However, it is possible to see by direct inspection 
for all the cases under consideration in the present paper, that those are indeed 
invariant w.r.t. the dynamics generated by the one-particle Hamiltonian. An¬ 
other nontrivial question arises in introducing the chemical potential g in order 
to investigate the thermodynamics of open systems. Our ansatz is to take for 
the occupation number n e at energy e and chemical potential /i (independent 
on the energy levels), 

1 

Ue ~ e /3(e)e—At _ 1 ' 

Another reasonable choice would be 

1 

77-6 e P(e-fi)(s-IJ.) — 1 

Both choices leads to equivalent results in the usual equilibrium situation (i.e. 
if j3 = const.), but these might lead to different situations in non equilibrium 
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ones. The most general situation in introducing the chemical potential would 
be that it also depends on the energy levels. Another quite natural choice will 
be briefly discussed in Section 7. 

As we are going to see, for the Bose case and also for g-particles with q > 0, 
we can exhibit a multitude of quasi-free states satisfying LEP and exhibiting 
BEC. Among such states, we will find many ones which exhibit unexpected prop¬ 
erties, even in the equilibrium situation. Due to Pauli exclusion principle, the 
Fermi case and also g-particles with q < 0, leads to only one stationary quasi-free 
state satisfying LEP. In fact, the same analysis can be done mutatis-mutandis for 
Fermi particle for which the condensation never occur. Consider the Canonical 
Anti-commutation Relation algebra which is the C* -algebra CAR(f}) generated 
by annihilators c(/) satisfying 

c{f)c\g) + c\g)c{f) = (f\g)t f,g£ (i. (3.8) 

For the gas of free Fermions, fix a density p > 0 and consider the (unique) 
solution of the equation 

/' d N(s) = 

J R+ + 1 P 

in the unknown g G R. Here, the cumulative function N is the so-called in¬ 
tegrated density of the states and describes the density of eigenvalues of the 
Hamiltonian in the infinite volume limit. At the knowledges of the authors, 
the integrated density of the states might not exists for a chosen finite volume 
exhaustion of the underlying physical space, or it may depend on the chosen 
exhaustion even in the amenable cases. The reader is referred to [17] for the 
rigorous definition of the integrated density of the states and some relevant prop¬ 
erties for some relevant amenable and non amenable cases. Then the quasi-free 
state determined by the two-point function 

^A^if)c(g)) = (g\(e^ h -» + 1)- 1 ]/), f,g g b 

is the unique quasi-free state fulfilling the LEP for the inverse local temperature 
function /3 and the chemical potential g. 

We end the present section by pointing out that our framework of local equi¬ 
librium can be easily generalised to self-adjoint, not necessarily semi-bounded 
Hamiltonians h, provided that for the local temperature function /?, f3(x)x is 
semi-bounded, almost surely w.r.t. the measure class determined by the resolu¬ 
tion of the identity e >->• eh{e) of h in (3.5). 

4 local equilibrium principle and the Bose-Einstein 
condensation 

In order to exhibit states describing the BEC even in the more general context 
of local equilibrium, we specialise the matter to the simplest model describing 
non relativistic free Bosons living on R d . To simplify, we put m = 1/2 for their 
mass. Analogous considerations can be done for Bosons on lattices Z d . Indeed, 
fix the functions in the class 


i)(R d ) C S(R d ) C L 2 (R d , d d x) 



made of the Fourier anti-transform of all the infinitely often differentiable func¬ 
tions with compact support in momentum space. The corresponding one- 
particle Hamiltonian will be 


j=i 3 


given by the opposite of the Laplace operator on R d . It is immediate to see 
that (i) before Definition 3.1 is satisfied. The one-particle Hamiltonian —A is 
nothing but the multiplication for the function 

fc 2 := ST k) 

3 = 1 


in the momentum space after Fourier Transform, where as usual, k = (k \,..., kd). 
Thus, we can consider any non negative Borel function {3 : (0, + 00 ) (0, + 00 ) 

with (3 > 0 almost everywhere w.r.t. the Lebesgue measure on R. 

The first interesting phenomenon is that, in the setting of local equilibrium, 
particles can condensate also on excited levels of the energies. The main result 
of the present section concerning the condensation regime for which p, = 0 (cf. 
Section 5) is summarised in the following 

Theorem 4.1. Let edlp2 ) P 2_ 1 € Ll oc (R d ), and the function x/3(x) € L^ C (R + ). 
Suppose that for some xo G [0, + 00 ), lim^^^p x/3(x) = 0. For each D > 0, con¬ 
sider the point-mass measure vdm '■= DS^. Then the quasi-free state u>p tVD k € 
S(CCR(T>(R d ))) with two-point function 


,f(p)g(p) d d p+Df{k)g(k) t f,g £ T>(R d ), (4.1) 


u^nAaHfWg)) := e/3(p2)p2 _ 1 

satisfies the LEP w.r.t. the local inverse temperature function /3 and chemical 
potential p = 0 , provided k 2 = Xq- 

Proof. Thanks to ed(p 2 j p2 —- G L 1 1 oc (R d ), (4.1) is well defined for each f,g G 
D(R d ). In addition, as xj3(x) G L^ C (R+) then 


0 2/3 (p 2 )p 2 


I/(P)| 2 


d d p < e 2 ll X / 3 ( X )r[°,r 2 ]ll 


I 2 

100 


cl d p 


/ p< r 


e 0(p 2 )p 2 _ l 


d p < +00 , 


Jr* e?(p 2 )p 2 - 1 

where supp(/) C D ri D r being the disk of radius r centred in the origin. Finally, 
if k 2 = Xq , the function e^ h l h f is uniquely defined in k as 

(e»/)(k) = ( lim e^ p2)p2 )f( k) = ( lim e 0(x)x )f( k) = /(k), 

' p—s-k ' v x—¥X0 ' 


because the function e^^ p ) p coincides a.e. with a measurable function which is 
continuous in k 2 = Xq. Collecting together, we have first that e^^ h f is in the 
domain of the form (4.1). In addition, by using the commutation relation (3.1), 
we compute 

w /VD,kK'7)a t (/)) = + e f3(p 2 ) P 2 _ ! ) /(P)g(P) d d P + D /( k )g( k ) 

0{p 2 )p 2 

IHp 2 )p 2 _ 1 /(P)g(P) dd P + De 0{k2)k2 f{k)g{k) = up >VD k (a f (f)a(e 0{h)h g )), 
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that is (3.7) is satisfied for p = 0. 


□ 


The two boundedness conditions in Theorem 4.1 have a different meaning. 
The second one x/3(x) € £&(»+). automatically verified in the usual equilib¬ 
rium setting when (3 = const., can be relaxed case-by-case when one consid¬ 
ers special examples of temperature function (3(e). The first one j £ 

LL0* d ) has an important physical meaning. First of all, we note that —^4—- £ 
oc m is equivalent to -^4—- £ 74(14) in the equilibrium situation. This is 
nothing but p c (/3) < +00 where p c ((3) is the critical density for the function (3, 
which does not depend on the energy levels in equilibrium thermodynamics. A 
key-role to study the appearance of the BEC is to look at the critical density 
p c ((3) at inverse temperature (3 = 1/T, see [13]. In our situation (including the 
equilibrium case), the critical density p c ((3) for the function (3 is given by 


Pc((3) 


d d p 

e /3 (p 2 )p 2 __ l ' 


(4.2) 


We will also show in Section 5 the natural role played by the critical density 
in studying the condensation effects. In the more general situation of LEP, 
IjfrTC e L loc( Rd ) is in general weaker than p c (/3) < + 00 . This means 
that we can easily exhibit quasi-free states describing condensation effects, for 
which p c ((3) = + 00 . A similar phenomenon can happen in studying BEC in 
equilibrium thermodynamics for inhomogeneous systems. We show that this 
last phenomenon is of different nature than the analogous one described in 
[18, 19]. In order to do that, we start by looking at the local density of particles 
p u (x) of a quasi-free state oj. It is given by 

a4 x ) = ui(a^ (5 x )a(6 x )), 


where S x is the Dirac distribution centred in x £ M , provided that the r.h.s. is 
meaningful, otherwise it is infinite. For the state in (4.1), we get for the local 
density of particles 


A^, k ( x ) 


IMp) 


e /3(p 2 )p 2 _ l 


d d p + D|4(0) 


From this simple calculation, we conclude that the D (x) is homogeneous and 
assume the form 


Pt*>l 3,„ D k ( x ) — Pc((3) + Pcond(^p,v D)k ) ) 

where 

D f ^ 

Pcond(^ D , k ) = aS 1 - VOl(A) J A l dx(0) l dd * = D - 

Thus, the local density of particles D (x) of the state wp,D is finite if and 
only if the critical density (4.2) of the model is finite too. This means nothing 
else that the models for which the critical density is infinite, still provide states 
exhibiting BEC effects, but all of them have infinite local density of particles. 
Namely, such states are in some sense unphysical. 

We pass to investigate the role played by the natural action an, R £ 0(d) 
of the rotation group on CCR(£>(R d )). For non trivial R £ 0(d) and non zero 
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k as in Theorem 4.1, we have i?k ^ k but |jl?kj| 2 = Xq. Thus, if there exists 
Xq > 0 such that lim x ._ >x . 0 /3(x)x = 0, the rotation symmetry is spontaneously 
broken. With v any bounded positive Radon measure on the sphere Sfc C 
of radius k, it is easy to see that for f,g € 2)(R d ), the quasi free-state with 
two-point function 


up, v ( a Hf)a(9)) ■■= J l ^]p 2 P ^ 1 d d P + /(p)s(p) d^(p) (4.3) 

satisfies the LEP. We immediately see that 

Pw 0 ,A x ) = Pc{/3) + v(R d ), 

where the last addendum on the r.h.s. takes into account of the portion of the 
condensate. In general, wp >v is not rotationally invariant. However, we get 

Proposition 4.2. With D > 0 and f lu the normalised rotationally invariant 
measure on the sphere Sk, the states in (4.3) are rotationally invariant. 

Proof. Let R £ O(d) be an orthogonal matrix with R l its transpose one. Define 
fn(x) := f(R t x). We easily get f R = (f) R . By taking into account that the 
function -^ 2 ^ 2 —j is rotationally invariant, we get after an elementary change 
of variable, 


ijJy,DQ. k {a R {a) {f)a{g))) = up :D n k W {(f) R )a((g) R )) 


/(i? t p)t/(i? t p ; d 


e ^(p 2 )p 2 — 1 


d p + D / f(Rp)g(R t p) dfifc(p) 


>Sk 


= L eWW-1 d " P + ° X /(P)5(P) d ^ (p) 

=up,DSi k {(J{f)a{g)). 


□ 


We already have discussed that the existence of quasi-free states exhibiting 
BEC for which the local density is finite, is determined by the convergence of 
the integral in (4.2) describing the critical density. To see that we can have 
the condensation even for spatial dimensions different from the usual one d > 3 
(or d > 2 for the fotonic/phononic Hamiltonian h( k) = At) , we consider simple 
examples for which (3 G C((0,+ 00 )), and 

In x 

f3(x) « x a ° for x —» 0 + ; f3{x) « a <*,- for x — > +00 . 

x 

The condition p c (/3) < +00 leads to > 0, whereas lim^o (3{x)x = 0 leads 
to «o + 1 > 0. Concerning the critical density, we compute for the one-particle 
Hamiltonian of the form h( k) = k s , s > 1 (to avoid unphysical models, even if 
the last restriction plays no technical role), 


r+00 


d d p 


e 0{p s )p‘ — l 


p d-l-«(a 0 + l) dp + 


r+00 


dp, 


11 




which converges if and only if 

s(a 0 + 1) < d < sctoo . 

This means that for the NESS considered in the present paper, it is possible to 
have or not BEC of free Bosons on for dimensions different from those one 
finds in equilibrium thermodynamics. 


5 a new approach to the condensation 

The standard way to find states exhibiting BEC on a locally compact manifold 
M, which can be (see e.g. [13]) or an infinitely extended network (cf. [18, 
19]), is to start from the Bose-Gibbs grand canonical ensemble (cf. [28]) of the 
finite volume theories based on a fixed exhaustion {A„} n6 N, A n f A/, together 
with the associated sequence of finite volume Hamiltonians {ffA„}neN- One 
considers the finite volume density 


PA(/?,Ai) 



c\N A (e) 

e 0(s)s _ 1 


associated to a compact region A C M. Here, with H A = f £de A (£), and N A is 
the cumulative function obtained as 


N A (£) 


Tr(E A ((—oo, e])) 
vol(A) 


which is well defined provided H A has compact resolvent. Under the usual 
conditions imposed to the function /?, the standard properties of the Hamiltonian 
H and the associated finite volume ones H A , the function p A (/3,/i) has the 
behaviour described in Fig. 1, with £o(A) > 0 the finite volume ground state 
energy of the system. This means that, for each fixed density p > 0, the equation 

f dN A (£) 

J R+ - l p 
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in the unknown p has a unique solution /z(A) G (—oo, /3(£o(A))£o(A)), which 
allows to determine the sequences |/i(A„)} of the finite volume chemical poten¬ 
tials. The critical density is now defined as 


Pc(P) ■= 


f d N(e) 

Jm + eM* - 1 ’ 


where N is the Integrated Density of the States associated to the infinite volume 
Hamiltonian of the system (see [17] for the rigorous definition), which exists for 
most of the models of interests with natural choices of exhaustions. If p c < +oo 
and the fixed density in Fig. 1 of the system p > p c , p(A n ) —»• 0 (possibly by 
passing to a converging subsequence if necessary), and the condensation phe¬ 
nomena take place. The standard approach briefly outlined here and commonly 
used to construct states exhibiting condensation, even if has a very clear phys¬ 
ical justification, is extremely sensitive of the boundary conditions imposed to 
the finite volume Hamiltonians (see e.g. [13] and the literature cited therein), 
and a-priori also of the chosen finite volume exhaustion {A„}„ e pj (see e.g. [38] 
and the literature cited therein). Due to the new phenomena relative to the 
BEC arising from the LEP described in Section 4, it is expected that this ap¬ 
proach might be not flexible enough to select/construct explicitly all the possible 
states exhibiting BEC. For example, if one fix the density as described above, 
one might select only rotationally invariant states which exhaust all of them in 
the equilibrium situation but not in the more general setting of local equilib¬ 
rium (cf. Theorem 4.1). In addition, even in equilibrium thermodynamics, we 
will find states exhibiting BEC but for which the density of particles is infinite 
(cf. Proposition 6.1). This suggests new approaches to the investigation of the 
BEC, which is the goal of the present section. We will see that this approach 
provides new and unexpected examples of states exhibiting condensation (cf. 
Section 6), even in the usual situation of equilibrium thermodynamics where 
the temperature is fixed. 

For our purposes, following the model described in Section 4 we consider 
h( k) = k s . In order to exclude unphysical models, we suppose s > 1, including 
the massive non relativistic Hamiltonian h( k) = k 2 (with the normalised particle 
mass m = 1/2), and the massless relativistic one h( k) = k (corresponding to 
the normalised speed of the light c = 1, or the velocity of the sound 1 in the 
case of phonons). In order to avoid technicalities, we make some reasonable 
restrictions to the function /3 : M + —► R + . Put fi(x) := fi{x)x and assume 

(i) ft G C([0, Too)), with infm + ft = 0. 

(ii) If E := /3 _1 ({0}), then E C [0, Too) is made of at most a finite numbers 
of points, with the empty set allowed. 

(iii) We assume also that 



i — i 

e 0(x) _ l 


da; < Too, 


for each a > ma xE. 

Assuming the continuity of temperature function is just to avoid technicalities, 
yet providing nontrivial examples, whereas infn + ft = 0 can be assumed without 
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loosing generality just by passing to the new function /3 — inf f3. The second 
condition is to avoid too general situations, perhaps already considered in Sec¬ 
tion 4. The third one says that e ^ (k s) k s_ 1 € L 1 1 oc (K d ) is simply equivalent to 
e T 1 (R d ) as for equilibrium thermodynamics. Denote 

S r := {(p, k) G R d x R d | p 2 = k 2 = r 2 } . 

the sphere of radius r living in the diagonal of the product space x R d . 
As our approach is quite general, we can manage the general g-Connnutation 
Relations for q £ [—1,1], ±1 being the Bose-Fermi alternative. By using Fourier 
Transform, we can consider the density creators and annihilators in momentum 
space which are operator-valued distributions, by putting 

a Hf)= [ /(k)a f (k)d d k, a(g)= [ g{k)a(k) d d k. (5.1) 

JR d JR d 

Then the ^-Commutation Relations, corresponding to (3.1), (3.8) for the Bose- 
Fermi cases, can be rewritten as 

a(k)a^(p) — qat(p)a(k) = 6(k — p)l. (5.2) 

It can be proven that the commutation relations (5.2), still generate an abstract 
C*-algebra, even for q £ (—1,1), see e.g. [24], In all the situations q £ [—1,1] 
considered here including the Bose/Fermi situation q = ±1, we denote by 
CCR 9 ((i) the involved C*-algebra. In addition, the quasi-free states are also 
well defined for the deformed cases — 1 < q < 1 because their 2n+ 1-point func¬ 
tions is 0, and 2n-point functions are described by the so-called q- determinant 
(being the case q = — 1 indeed a determinant, known as the Slater determinant) 
which is seen to be positive for the case considered here, see [11, 31, 34]. The 
LEP for g-particles, q £ [—1,1], still assumes the form (3.7) in Definition 3.1. 

The new approach in searching quasi-free states u> exhibiting condensation is 
to look at those for which the two-point function is given by a distribution F, b, 
at least for the situation of free particles considered in the present paper. For 
such a purpose, we use the conventional integration symbology of the theory of 
Distributions (cf. [40]) for the natural pairing F(f) between smooth functions 
/ and elements F of the topological duals 

/ x F £ R(K Z ) x D'(R‘) F(f) =: [ F(|)/(|) d'C • 

Jr 1 


With this symbology, for f,g£ 2)(K d ), we compute by considering (5.1), 


uj(a\f)a(g)) = f w(a t (p)a(k))/(p) 3 (k) d d pd d k 
JR d xR d 


/R d xR d 


F u (p, k)/(p)<?(k) d Pd k, 


where F u £ D'(]R d xK d ) is some distribution given by F u ; (p, k) := w(a^(p)a(k)). 
By positivity, such a distribution should satisfy for each / £ D(R d ), 


F UJ (p 1 k)f(p)f(k)d pd k > 0. 
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A positive definite distribution as above is said to be a kernel. In addition, a 
kernel is automatically real, that is it should satisfy in the sense of distribution, 


k,p) = P,k), 


(5.3) 


where the bar stands for complex conjugation. 

Let (jJ £ §(CCR 9 (i)(K d )) be a quasi-free state whose two-point function 
satisfies (3.7). Suppose further that its two-point function is given by a kernel 
F u as described above. When we try to impose condition (3.7) to F u , a product 
of a distribution with a function which is in general not smooth, shall appear. 
This product does not define any distribution in general. Thus, we suppose 
further that F u ( p, k) defines still a distribution. Notice that, by 

positivity, e^^ k ^ h ^~^F u (p,k) also defines a distribution. This will play a 
crucial role in the following. By using the commutation rule 5.2, the LEP (3.7) 
leads for the chemical potential /i, and the parameter q £ [—1,1] to 


/R d xR d 


(e /3(M k ))M k )-M )a;(o t (p)a(k))/(p)g(k) d d p d d k 


=w(a t (/)a( 7 /i (/i)^)) = w(a(g)a T (/)) = (g\f) + gw(a T (/)a(g)) 


(<5(p - k) + <7w(a t (p)a(k)))/(p) 5 (k) d d p cl d k, 


obtaining 


[ L/J(h( k ))M k )-P _ q) Fu ( p, k) - 6 (p - k) 

J R d xR d L 


/(p)ff(k) dp d d k = 0 . 


As the last should be satisfied for the total set in CD(IR d x R. d ) generated by 
elementary tensors, we obtain 

( e /j(h(k))h( k )-,. _ q )F ul (p, k) = S(p - k). (5.4) 

Combining the reality condition (5.3) with (5.4), we get also 

( e /J(/*(p))fc(p)-»* _ q ) Fuj ( p, k) = 6(p - k). (5.5) 

The above computations tells us nothing but that the kernel F^ reproducing a 
quasi-free state on CCR g (i)(K d )), which solves (5.4) (or equivalently (5.5)) are 
the natural candidate to describe those satisfying the LEP for the g-deformed 
relations at the inverse temperature function 13 and chemical potential g. As we 
will see below, in searching solutions of (5.4), we should divide for the function 
e (/3(/i(p))ft,(p)—/i) _ q w ] 1 i c ] 1 by positivity should be everywhere positive, with zero 
possibly allowed on a negligible set w.r.t. Lebesgue measure on the momentum 
space. This allows to compute the possible range of the chemical potential, 
leading to g G (—oo ,g q \, where fi q = +oo for q G [—1,0], and g q = — lng for 
q G (0,1], 0 being the limiting value between the two situations. As it will be 
more clear below, the the previous cases can be considered as the Fermi/Bose- 
like alternative, with the separation case q = 0 is known to be the the Boltzmann 
(or free) one. 

In the Fermi-like/Boltzmann cases q £ [—1,0], including the limiting Boltz¬ 
mann case q = 0, we can freely solve (5.4) for each g £ R. So the concept of 


15 



critical density plays a role only in the Bose-like case q £ (0,1]. For such values 
of q, (5.4) (or equivalently (5.5)) can be solved only ii p < /. i q := — lng, being 
/ i q the critical value of the chemical potential for which the condensation can 
occur. Concerning the critical density at the inverse temperature function /?, it 
is defined as p { 'f\f3) := p^(/3,/r 9 ) which leads to 



d rf p 

P(h(jp))h(p)-it q 


- q 


P { c\P) 

q 


(5.6) 


We will se below the role played by the limiting value of the chemical potential 
p q and critical density pi?\ in the appearance of the condensate of g-particles. 
As before, pc is simply denoted as p c . We also omit the dependence on /3 when 
this causes no matter of confusion. 

According to the previous computations, the following theorem explains the 
form of the kernels which are the candidate to reproduce the two-point function 
of a quasi-free state satisfying LEP. Among the other things, we show how the 
condensation regime naturally emerge without using the thermodynamic limit 
of finite volume theories. 


Theorem 5.1. Suppose that the inverse temperature function (3 fulfils (i)-(iii) 
above. Let F be a kernel on such that e W h ( k )) h ( k ) _#1 )i i ’( p, k) is still a dis¬ 
tribution satisfying (5.5). Then the following assertions hold true. 

(i) If p > — ln(0 V q) (with the convention that — InO = +oo,), no of such 
kernels can exist. 

(ii) non condensation regime: For each p < — ln(0 V q), there exists only 
one kernel as above having the form 

F(k, p) = e/S(h J ))h(p) _l _ q • (5-7) 

(Hi) condensation regime; Let q £ (0,1] and p = p q . If E = 0, then 
pi q ' > < +00 and 

k ) = q ( e PWp))h(p) - 1) ( 5 ' 8 ) 

is the unique kernel fulfilling the hypotheses. 

Let E ^ 0, and suppose that /3[e. + \e is infinitely often differentiable. 

If pi 9 \/3) = + 00 , there is no of such kernels satisfying the additional 
condition |/| < |.g| => F(f <g> /) < F(g <g) g). 

If pc < + 00 , then F assumes the form 

F ( p ’ k ) = g ( e /3(h(p))h(p) _ x) + G ( p ’ k ) ’ ( 5 ‘ 9 ) 

where G £ D'(R d x K d ) is supported in (J {§ 1 | r £ E}. 

Proof, (i) and (ii): Being n p = e ft(h( P ))h(r,)-^_ q the density of occupation number 
at momentum p and chemical potential p, it must be a.e. positive, w.r.t. the 
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Lebesgue measure. This excludes p > — ln(0 Vg). Conversely, if p < — ln(OVq') 
we can uniquely solve (5.4) by obtaining (5.7). 

(iii): By (ii), the condensation regime can occur only if q £ (0,1] and p = 
— hi q. So we reduce the matter to this case. If E = 0, then condition (iii) 
at the beginning of the present section together (5.6) imply that p/P < +oo. 
In addition, being e^Pd p ))fr(p) — 1 > 0 everywhere, (5.5) can be freely solved 
w.r.t. F giving as the unique solution (5.8). Suppose now E ^ 0. We first show 
that, under the quite natural conditions imposed to /? and to the kernel F, 
the critical density must be finite. Thanks to (5.6), we can reduce the matter 
to q = 1. Suppose that p c = +oo. By condition (iii) above, it means that 

Jb pwx-i d® = +oo, at least for some small neighbourhood B Xo centred in 
Xo £ E of radius 5 containing eventually no further points of E. Choose a 
positive function ip £ D(R d ) which is identically 1 on the spherical shell of 
thickness (5 around the sphere h{ p) = xq (the sphere of radius S/2 centred in 0 
if Xq = 0), and whose support does not contain further points p with h( p) € E. 
Choose a sequence C ©(M d ) of positive functions with p n < ip, such 

that the sphere h{ p) = Xq is not contained in their support, monotonically 
converging point-wise a.e. (w.r.t. the Lebesgue measure) to (p. By hypothesis 
on F, we have for each n £ N, 

f +(p,k)</? n (p)v?„(k)d d pd <i k < / F(p. k)^?(p)<p(k) d d p d d k < +oo . 

Jm. d xR d J R d xR d 

On the other hand, as Jit/PPp/h/lPi is a smooth compactly supported function, 
by using (5.5) and the Monotone Convergence Theorem, we obtain 


/ R d xR d 


/ R d xR d 


F{ p, k)<p„(p)<p n (k) d d p d d k 

' p 0(h{p))h{p) _ i \ pf U'A +”(P)+"(^ : ) id idl 
e iKlP,^ ej8(Mp))h(p) _ 1 d Pd i 

<5(p-k) 


Jwi d xM. d eW(p)Wp) - l 

/' Vn( P) 2 


> 


e/5( ft (p)) ft (p) — l 

f d d p 

j B e P(h(p))h(p) _ ^ 


d d p 


+n(p)+n(k) d d pd d k 

/• ^(p) 2 


J R d e /3(h(p))h(p) - 1 

d d p = +oo 


dp 


which is a contradiction. 

Suppose now that p c < +oo. Then defines a distribution on 

x because of (iii) at the beginning of the section. Then 


G(p,k) 


F( p,k) 


<5(p - k) 

e 0(h(p))h(p) _ l 


defines also a distribution. Now we show that supp G = U {§ 1 | r £ E). Pick 
ip £ ©(K d x K d ) with K := slippy C {§ i | r £ E} c . We can suppose without 
loosing generality that K is a smooth manifold with boundary. Consider the 
compact slices Ki := K D {(p.k) | h{ k) £ E}, K 2 := K n {(p, k) | h( p) £ F} 
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and choose open neighbourhoods V D K, , i = 1,2, such that the distances 
between K\(\\ U Vf) and A^, i = 1,2, AT\(Vi U V 2 ) and (J {§^i | r G A} is 
greater of 2d for some S > 0. Then there exists a finite open covering of I\ made 
of Vi D Ki, i = 1,2, (which cover the AT.;), and some balls centred in 

some points {£j}” =1 C K\(V\ U V 2 ) of radius d (which cover K\(V 1 U I4)). Fix 
a smooth partition of the unity 

n 

ip = V >1 + 1 P 2 + ^2 fi 

i=l 


subordinate to this covering. After using (5.5) for the addendum containing 
ipi, (5.4) for that containing ip 2 , and finally indifferently one of them for the 
addenda containing the fi, we first note that 


l(Pik) 

e 0(h(p))h(p) _ l 


Vh(p,k) + ^Tfi(p,k) 


tp(p 1 k)V> 2 (p 1 k) 

e /3(?i(k))/i(k) _ l 


are smooth compactly supported functions. Then we compute, 


G(p,kMp,k) d d p d d k 


(e 


0(h(p))h(p) 


- l)-F(P,k) - d(p - k) 


x e /3(X)WP) - 1 (^ l(p ’ k) + E M P> k )) d<i P d<ik 


( e ^(h(k))/,(k) _ i)F(p, k) - d(p - k) 

R d xR <i 

^(p,k)y> 2 (p.k) d d _ 

x e /3(h(k))/ l (k) _ x d Pd k °- 


□ 


We point out that the investigation of the condensation phenomena of q- 
particles in the setting considered in the present section, can be reduced to find 
the positive definite solutions of the equation (5.4) (or equivalently of the twin 
one (5.5)) in the space of the distributions. As shown in Theorem 5.1, it is a 
difficult task to find all the possible solutions of that equation. In addition, not 
all the distributions G appearing in (5.9) give rise to a quasi-free states satisfying 
LEP and exhibiting condensation. However, we will see in Section 6 that, even in 
the standard equilibrium case when f3 is the inverse temperature of the system, 
(5.4) admits solutions of the form (5.9) with nontrivial G, hence describing 
condensation effects, which are completely new and unexpected. Concerning 
the quasi-free states described in the previous section, the two-point function 
(4.3) also have the form (5.9) with nontrivial G, including those in Proposition 
4.2 which are rotationally invariant, even if the local inverse temperature is more 
general than those considered in Theorem 5.1. 


6 new states describing condensation 

In the present section we show that new states exhibiting condensation can occur 
even in the standard equilibrium thermodynamics. We reduce the matter to the 
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Bosonic case with /3 the fixed inverse temperature. As usual, we put h(p) = p s , 
s > 1 , and reduce the analysis to the simplest case of the partial derivative along 
only one direction k\ in the momentum space. As for ip £ T>(R d ), e‘ p *<p(p) is 
in general not smooth when s = 1 or s is not integer, we slightly extend the 
space of the test functions, even if it is not necessary for the unique physical 
case s = 2 for which it is possible to exhibit such new states describing BEC. 
We put 

f)o = span {e lht (p \ ip € D(K d ), t £ R} . 

In this case, the one-particle dynamics is meaningful on f)o by construction. 


Proposition 6.1. Fix D > 0 and suppose that the critical density p c < +oo. 
Then the states whose two-point function is given by 

"<“ t(/ >< s >> = L + + < 6 - 1 > 

are equilibrium states for the dynamics generated by the one-particle Hamilto¬ 
nian h( p) = p s , provided s > 1 . 

Proof. Fix a generator / = e iht <p. It is straightforward to see that 


f | e ^ 7( P )| 2 

eP pS — 1 


d d p < +oo. 


In addition, if s > 1 we get 

^ (P) = /3sp s - 2 Pl e^ 3 f(p) + e^f^(p). 
dp i dpi 

Collecting together with y(p) = e^ p °, we conclude first that w(a t ( 7 /)a( 7 /)) < 
+oo, which means that 7 / is in the domain of the form describing the two-point 
function of to. Second, d{ Q p ^ (0) = J^-(O), leading to (3.7). Namely, the states 
in ( 6 . 1 ) are equilibrium states at inverse temperature /? and chemical potential 

p = 0 . □ 

The two-point function of the usual state p exhibiting the condensation of 
free Bosons is given by 

p(a\f)a(g))= [ d d p + Df(O)^O ), (6.2) 

J Rrf ePP - 1 

and has the form (5.9) with 

G v (p,k) = <5(p)<5(p — k). 

Despite the product of two delta distributions, G v { p, k) makes sense as a distri¬ 
bution, giving precisely the condensation portion after the application of Fubini 
rule in the formal integration as explained below. Concerning the two-point 
function in (6.1), it assumes the form (5.9) with 

r f , X _ d' 2 (6(p)6{p-k)) _ ff 2 G v { p,k) 

’ dpidhi dpidki ’ 
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which leads to 


G u { p,k)/(p)5(k) d p d k 


d 2 G v { p,k) 


T <p 

dpidki 

r dl ,%( k ) 


/( P ) 5 (k)d d Pd d k 
(5(p - k) 




The higher (even) derivatives of the delta-distribution can give rise to a con¬ 
densate distribution only if the power s in h( p) = p s is sufficiently big. Con¬ 
versely, it never appear in the case of photon/phonon one-particle Hamiltonian 
h( p) = p, and only states like (6.1) involving the second derivative of the delta- 
distribution and not the higher ones can appear in the free massive Bosons for 
which h( p) = p 2 . 

The local density of particles p w (x) of the states in (6.1) is finite everywhere 
in the configuration space: 


Pw(x) := w(a* (6 x )a(5 x )) 


d rf p 

e 0P s - 1 


+ D 



+ Dx\. 


(6-3) 

Conversely, the mean density is infinite due to the contribution of the conden¬ 
sation term -D|||y(0)| 2 = Dx 2 . As the standard procedure in constructing 
states exhibiting BEC as infinite volume limits of finite volume theory is to fix 
the mean density of the model, the last computation easily explains because 
such terms connected to the second derivative cannot appear in the standard 
investigation based on Bose-Gibbs prescription for the grand canonical ensemble 
construction. 

States which are rotationally invariant (i.e. isotropic) are easily given by 


v( a Hf)a{g))= f d d p + -PVp/(0)-V p g(0), (6.4) 

as well as states exhibiting condensation, which are connected to higher deriva¬ 
tive of the Dirac distribution can be easily wrote down, provided that s is suffi¬ 
ciently big. It is also possible to exhibit such states involving higher derivatives 
of the delta-distribution in local equilibrium, obtaining a richer situation. 

For all the states described in (6.4) (and also for those in (6.1)), the action of 
the spatial translations in the configuration space T x is spontaneously broken. 
We easily compute 


w o T x (a f (/)o(ff)) = u(T x ((J(f)a(g))) 

= u(a\f)a(g))+D[x 2 f{0)g(0 ) + *x-( 5 (0)V p /(0) - /(0)V p £f(0))] , 

where as usual the dot denotes the real inner product. Contrarily to the states 
in (4.1) for the rotation symmetry, for the states in (6.4) also the local density 
of the particle changes for a term relative to the condensate which is simply 
computed as 

PwoT y (x) = Puj (x) + Dy 2 . 
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Figure 2: The density p^ cond ) ^ 2 ) in momentum space relative to the states 
(6.2) for two values E\ > £ 2 . 


Fix any mollifier <5 e (p), converging in the sense of distribution to the delta 
distribution <5(p) as e 4- 0. In order to have an idea of the distribution of the 
condensate in momentum space of the usual situation (6.2) and that of (6.4), 
we end the present section by drawing the projection in the fci-fc 2 plane of 

p^ ond) (k 1 ,k 2 )=S e (k 1 ,k2,0) 2 , 

P ^ d \k ly k 2 ) = i|v p 4(p - k)r P=0 , fe3 =o|| 2 , 

where ui is the standard equilibrium state exhibiting condensation (6.2) (cf. 



Figure 3: The density pu° nd \ki, ft 2 ) in momentum space relative to the states 
(6.4) for two values £1 > £ 2 . 

Fig. 2), and one in (6.4) (cf. Fig. 3), respectively. 


7 remarks and outlook 

In order to study the condensation regime, we have met the problem concern¬ 
ing the way to introduce the chemical potential. We briefly discuss a further 
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reasonable choice 


n, 


(7.1) 


1 


1 


e Z~ 1 eP( e ') £ — 1 e /3(£)(e-/i) _ X 

corresponding to choice for the the activity 

Z {/3, p) = ef** = ePW* . 


The analogous equation to (5.5) assumes the form 

( e 0(Mp))(Mp)-#O _ q )F u (p, k) = J(p - k). (7.2) 

It is simply to see that this different possibility introduces no change for the 
cases q £ [—1,0], and for the limiting Bosonic case q = 1. For the Boson-like 
cases q £ (0,1), we compute for the sake of completeness, the critical value 
of the chemical potential p q , and hence the relative critical density pfP • The 
key-point is the following 

Lemma 7.1. Let the positive function fd : (0,+ 00 ) —> [0,+ 00 ) be continuous 
and vanishing of a set E C (0, + 00 ) such that the set L (E) of its cluster points 
is discrete, //limsup^m < + 00 , then the function G : R+ —» K+ given by 
G(p) := sup (M {/3(x)(p - x)) is strictly increasing and continuous. 

Proof. Define /3(0) := limsup^.Q and fix p > A. There exists a x £ [0, A] for 
which /3(x) > 0 such that by applying Weierstrass Theorem, we get 

G(p) = sup (f3(x)(p — x)) > sup (/3(x)(/z — x)) 

(0,/d (0, A] 

> sup (/?(x)(A — a:)) + (3(x)(p — A) > sup (/3(x)(A — a;)) = G( A). 

(0,A] (0, A] 


Thus G{p) is strictly increasing. Concerning the continuity, first notice that 

lim sup (f3(x)(p — a;)) = sup (/?(x)(A — a:)), 

/4A (o,A] (0,A] 


because for x < A, 


b{x){p - x) = /?(x)(A - x) + (i{x){p - A). 


In addition, 

lim sup (f3(x)(p — x)) = 0 . 

*4 A [A,m] 

Collecting together, we get 

limG(/z) = lim ( sup (/3(x)(p — x)) V sup (/3(x)(// — x))) = G(A) 
m4- a a 4A ( 0jA ] [A,/*] 

For the reverse limit, first notice that, with p fixed, and A < p 1 

A i->- sup (/ 3(x)(p — x)) 

(0,A] 


is continuous. In addition, 


lim sup (j3(x)(p — A)) = 0 . 
At>(o,A] 
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Then we get with p > 0 fixed, A < p, and x G (0, A], 


/3(x)(p — x) — sup (Ji{x)(p — A)) < f3(x)( A — a:) < /3(x)(p — x ), 

(0, A] 


which leads to lim^ G{ A) = G(/z). 


□ 


The choice (7.2) for the occupation numbers introduces technical troubles 
only for the unphysical cases q G (0,1), for which we are going to determine the 
critical chemical potentials responsible of the condensation effects. 

Proposition 7.2. Let /? satisfy the hypotheses of Lemma 7.1. 

(i) Suppose that limsup x i 0 /3(x) = +oo, and define p q = 0. 

(ii) Suppose that limsup x , 0 0(x) < +oo. Then for p > 0, the function 


F q (p)'-= - 77 ^- 96 ( 0 , 1 ), 

SU P(0, M ] {p{ x Kt-x)) 


(7.3) 


is continuous and strictly decreasing. Define p q = 0 if the equation 

FM = 1, (7.4) 


has no solution which happens if and only i/Timsup^m @(x) < 1, or p q as 
the unique solution of (7.4). 

In both situations, the following hold true. If p < p q then 
( e /3(x)(x-^ _ g ) > s > o > xgK+j 


and if p > p q there exists an open interval I C R+ such that 

(g/JCrXx-M) _ q) < -s C 0 , x€l. 

Proof. If limsup^m /3(x) < +oo, Lemma 7.1 assures that (7.3) is a continuous 
strictly decreasing function. In all the situations, if x > p, then > 

1 > q. If p < p q and x < p then for some £ > 0, 

j3(x)(p — x) < /3(x)(p q — x) — e < sup (/3(x)(p — x)) — e =—lnq — e, 

(o ,Uq] 

which leads to 

e 0(x)(x-/i) > qe c >q . 

Conversely, if p > p q , we get in all the situations that there exists an open 
interval I G (0, p) such that for some £ > 0, 


/3(x) (p — x) > — In q + £, x G I. 


But this immediately implies for x G I, 

e P(x)(x-n) < qe -e K q 


□ 
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Notice that, being In q — 0, (7.4) is meaningful even for q = 1 providing 
pi = 0 as expected. In addition, if /3 is constant, then /3p q = —Ing, which 
reduces to the usual formula of equilibrium thermodynamics. Proposition 7.2 
explains why p q solving (7.4) is the critical value of the chemical potential for the 
g-deformed situation, q G (0,1), for the choice (7.1) concerning the introduction 
of the chemical potential. The corresponding critical density is easily given by 


P 


(?) = 


d d p 

e P(h(p))(h(p)-Vq) — q 


Unfortunately, there is no direct relation between { p ^ j q G (0,1)} and p^K 
This introduces some additional technical troubles in order to manage (7.2) for 
the unphysical cases 0 < q < 1 in the condensation regime p = p q , without 
affecting the substance of the results. We do not pursue this direction, which 
perhaps can be managed in details for particular choices of the local temperature 
function f3. 

We have discussed new ideas of Local Equilibrium Principle which is a suit¬ 
able candidate to select stationary states in non equilibrium thermodynamics. 
Roughly speaking, the Local Equilibrium Principle asserts that the tempera¬ 
ture is a function of the energy levels. Even if the Local Equilibrium Principle 
cannot be directly extended to general dynamical systems describing models 
with infinitely many degrees of freedom, as well as general condensed matter 
states appearing in nature, we have provided a wide class of nontrivial new 
and unexpected examples involving Boson particles and describing the Bose- 
Einstein Condensation. It is expected that the standard construction based on 
the Bose-Gibbs grand canonical ensemble involving the thermodynamical limit 
of the finite volume theories cannot directly cover all the examples we have 
found. This suggested a new approach based on the theory of the Distributions. 
This new approach allowed us to manage in an unified way the g-Commutation 
Relations, and the condensation of the so called g-particles can appear in a nat¬ 
ural way only for the Bose-like situation g G (0,1]. By using the new approach 
to BEC described in Section 5, we found new examples even in the standard 
case of equilibrium thermodynamics when the temperature is kept fixed by an 
external thermal bath. All these new states might find applications in high 
density/temperature non equilibrium physics, quantum optics, theory of super¬ 
conductivity. Among such possible applications, we briefly outline the following 
ones as promising perspectives. 

The first one concerns the states in the usual case of equilibrium thermody¬ 
namics described in Section 6. One consider the axially symmetric states whose 
two-point function is given by 


w(a f (/)a(g)) 


/(P)g(p) 

g(3p 2 /2m* _ ^ 


d 3 p + D 




concerning the distribution of some quasi-particle of effective mass in* at inverse 
temperature /? > 0. The spatial local density is easily computed as in (6.3) 


Pu,{x,y,z) 


L e /V/ln P ._i + D (x 2 + y 2 ). 


(7.5) 


The last addendum 

p c ° nd (x,y,z)=D(x 2 + y 2 ) 


24 



describes the spatial portion of the condensate, and has the form of a z-axial 
symmetric paraboloid. The helium superfluidity was explained (cf. [27]) by 
the excitations of the quasi-particles corresponding to the rotonic part of the 
spectrum of He 4 as briefly outlined in the introduction. Even if the superfluidity 
seems not directly connected with the BEC, it should be noted the surprising 
analogy between the spatial distribution of the condensate (7.5) and the profile 
(i.e. the meniscus) of the rotating He 4 superfluid, apart from the finer regular 
structure made of vortices, appearing inside the rotating superfluid, see e.g. 
[16]. 

The second one concerns the cosmological application of the Local Equilib¬ 
rium Principle for highly concentrated systems. For example, due to the net 
energy flow, the black hole evaporation (cf. [ 8 , 9, 23, 39]) is a non equilibrium 
process. Then we can find states which are far from the equilibrium, for which 
the temperature is a function of the total energy of the black hole which of course 
changes during the evaporation process. Another possible application might be 
to the thermodynamics of neutron, and also of exotic Boson and quark stars. 
For such systems subjected to extreme conditions, first it might be expected 
that NESS are the natural candidates to describe their thermodynamics. Sec¬ 
ond, due to high pressure, the Fermions might form BCS pairs like in He 3 or in 
superconductors, and the BEC might take place. On the other hand, the Local 
Equilibrium Principle allows states admitting portions of the condensate even 
in excited levels. In fact, states similar to those in (4.3) with 

h( p) = ca/p 2 + m 2 c 2 — me 2 

the one-particle Hamiltonian of a free Boson gas in a neutron or exotic star 
(with m being the mass of the involved Boson and c the speed of light, where 
in a rough approach we are neglecting the gravitational interaction between 
particles), describe a portion of the condensate in excited levels which cannot 
appear in the usual equilibrium thermodynamics. Such possible condensation 
effects in excited levels might partially explain the open problem of the dark 
matter of the universe. 

We conclude by pointing out that the last two possible applications of the 
Local Equilibrium Principle are only ideas which deserve of further insights, and 
are still very far to be understood at this stage. 
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